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ABSTRACT 

We account for all the image distortions relevant to weak gravitational lensing to second 
order. Besides the familiar shear, convergence, rotation and flexions, we find a new image 
distortion with two distinct descriptions, the twist and the turn. Like rotation, this distortion 
is not activated gravitationally to first order, but will be activated by systematic effects. We 
examine the rotational properties of twist and turn, and their effect on images in real and 
shapelet space. We construct estimators for the new distortion, taking into account the centroid 
shift which it generates. We then use these estimators to make first constraints on twist using 
the STAGES HST survey; we find that the mean twist estimator is consistent with zero. We 
measure correlation functions for our twist estimator on the survey, again finding no evidence 
of systematic effects. 

Key words: cosmology: observations - gravitational lensing. 



1 INTRODUCTION 

Weak gravitational flexion is a relatively new addition to the 
panoply of gravitational lensing effects, but has considerable 
potential for measuring substructur e in the density distribution 
of matter in the Univers e (see e .g. Goldberg & Nataraian 2002^; 
Goldberg & BaconI l2005l: llrwin & Shmakova 2005 ; Bacon et al.i 



20061 : lOkura et al. 1120081 : 1 Schneider & Eill2008h . 



Flexion is proportional to third angular derivatives of the pro- 
jected gravitational potential along the line of sight. As such, it 
is at the next order of differentiation compared to shear and con- 
vergence, which are the more s tudied weak lensing measures (see 
iBartelmann & SchneideJl200lL for an extensive review). Since, 
as we shall see, there are two independent combinations of third 
derivatives, there are two different flexion effects: the 1 -flexion, 
which is a vector distortion leading to objects being skewed; and 
the 3 -flexion, which is a spin three distortion changing circular ob- 
jects into trefoils. 

Up until now, these have been the only known image distor- 
tions at this order. However, in this paper we will show that there 
is a further neglected image distortion at the flexion level, with 
two alternative descriptions which we call twist and turn for rea- 
sons which will become obvious. This distortion is not activated by 
gravity under the most straightforward approximations; but it will 
be activated by systematic effects. The latter are of great concern 
to weak lensing, so finding a further signature of systematics is po- 
tentially very valuable to upcoming lensing surveys. 

In this paper we show how twist or turn arises, and account for 
why it has not been noticed before. We show how it affects images 
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in real and shapelet space, and give details of how it can be mea- 
sured with fairly straightforward estimators. We then measure twist 
for the first time using t he Space Telescop e A90 1/902 Galaxy Evo- 
lution Survey (STAGES. lGrav et al.l2007h . a large mosaic observed 
with the Hubble Space Telescope (HST). We show that the twist is 
consistent with zero for STAGES on all scales, both in terms of 
its mean values and its correlation functions, incrementally adding 
confidence in the management of systematics for this survey. 

The paper is organised as follows. In Section [21 we recount 
the theory of image distortions in weak lensing at the more studied 
first order. We note that there is already a non-gravitational mode 
at this order; image rotation. We write the distortions in terms of 
Pauli matrices, which will give us the necessary clues for how to 
treat higher order distortions later. 

In Section [S] we extend the account to second order. We find 
that there are combinations of Pauli matrices orthogonal to those 
describing the conventional flexion degrees of freedom; these or- 
thogonal combinations give twist and turn distortions. We are there- 
fore able to write down for the first time the complete weak image 
distortion to second order, and show how twist and turn are related 
to one another observationally. 

Section|4]describes the behaviour of twist/turn. The rotational 
properties of the distortion is worked out, and we find that is is 
a vector quantity. We show the impact of twist and turn on sim- 
ple images; we find that they do not affect the shape of circularly 
symmetric images, but only images with non-zero ellipticity We 
show explicitly the nature of twist and turn in shapelet space, prov- 
ing that they have no impact on circularly symmetric sources, and 
derive how they move power between shapelet coefficients. 

In Section[5]we go about finding practical estimators for mea- 
suring twist and turn. We derive simple estimators in shapelet 
space. Noting that like flexion, twist and turn affect the centroids 
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of objects, we correct the estimators by constructing slightly more 
complicated expressions which take this shift into account. How- 
ever, we will show that our estimators are not perfect; they would 
respond to flexion if it is present, and should therefore be treated as 
estimators of any second-order systematics, or of (real or system- 
atic) twist on scales where flexion is negligible. 

In Section [6l we use these estimators to constrain twist obser- 
vationally for the first time, using the STAGES HST survey. We find 
that our twist estimator has a larger variance than flexion estima- 
tors, and that its mean value is consistent with zero in the STAGES 
data. We measure correlation functions for twist estimators, again 
finding that they are consistent with zero systematic in STAGES. 
We summarise our results and conclude in Section[7] 



2 IMAGE DISTORTIONS TO FIRST ORDER 

We begin by discussing i niage distortions in weak lensin g to first 
order (for more details, see lBartelmann & SchneideJ200lh . We can 
describe the effect of lensing as a mapping between the surface 
brightness fs of a galaxy at a position (/3i , /32) in the source plane, 
and the surface brightness // at a position (^1,^2) in the image 
plane: 



fI{0^)=fs{f3^)=fs{A.J0J) 



(1) 



where we have set the origin of Oi and ft to the centre of light in 
the respective planes. A is the Jacobian matrix which maps image 
positions to source positions. 



dOi 



(2) 



For lensing with a single lens plane, and assuming the Born ap- 
proximation, this is given by 



■ didj^p 



However, our study of A is not complete. We have specified 
three quantities in A, i.e. k, 71, and 72. But A is a four element ob- 
ject, so there is a further degree of freedom which we have missed. 
We quickly realise that this is a rotation p, i.e. 
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(9) 



for small rotation angles p. Whereas k and 7 can be written as sec- 
ond derivatives of the lensing potential, this is not possible for p. It 
is not activated by gravity in our approximation (due to the inter- 
changability of the second derivatives of the gravitational potential, 
didjip — djdiip), but may be present in a real lensing survey as a 
sy stematic (see e.g. th e rotation caused by the telescope constrained 
bv Bacon et al.l200Qh. This rot ation has been described previously; 
seee.g. lHirata&SeliakI (l2003h . 

We will find it convenient to write A as a sum of Pauli matri- 
ces, as these provide an orthogonal basis for studying further de- 
grees of freedom at the next ord er of weak lensing appr oximation. 
The Pauli matrices are given by (lArfken & WebeduOOSh 
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SO we can write A as 



A = {1 — K,)I — 7icr3 — 72cri + pia2. 



(10) 



(11) 



(3) 



We will need one further concept: in a weak lensing context, it is 
usual to assume that the shear and convergence are small and con- 
stant across an object. We can then write the surface brightness 
mapping as 

fI{0^) = fs{S^J0J^{A,J-6^J)0J) c^ fs{0^)+{A,J-6^J)0J^^fs{0^Xl2) 



where ip is the lensing potential, i.e. the gravitational potential suit- 
ably projected into 2D. We can therefore write A as 



We will now modify this to show how flexion and the new distor- 
tions enter. 



A 
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with the convergence k, given by 

i^= 2(^1 +^2 )'0 

and the shear 7^ given by 



71 



>' 



■ 9l)^, 



72 = did2ip- 



(4) 



(5) 



(6) 



There is an alternative notation that is useful to us, introduced by 
iBacon et al.l (|2006). We define the complex derivative ^ = ^1 + 
id2 ; in cylindrical coordinates this is given by 



icf) 



d id 



with radial coordinate and azimuthal coordinate ( 

7 = 71 + 272, and then 



1 



ar^, 



1 



7 ^ 



ddil^. 



(7) 
. We also define 

(8) 



Besides simplifying notation, this format elucidates the spins of the 
quantities; when d is applied, the &^^ term in equation raises the 
spin by one. Similarly, the application of d* lowers the spin by one. 
So since '0 is a scalar, so is n, while 7 is spin 2. 



3 IMAGE DISTORTIONS TO SECOND ORDER 

The further step taken by flexion studies is to note that in reality, 
shear will vary across an object. If we keep A as a constant across 
the object, we need a further ter m in a Taylor expansion in the sur- 
face brightness map, as given bv lGoldberg & BaconI (l2005h . 



// 



fs { AijOj + -DijkOjOk 



(13) 



This introduces the D tensor; if we suppose that its components 
are purely due to a variation of A across the image, we can write 
Dijk = dkAij. Then by differentiating equation ilTl we find 



Diji -- 

Dij2 



-diKl - dijia3 - 

-82^,1 — ^271 CTs 



^172C^1 + dlpi(T2, 

- 82^20' I + d2pio-2. 



(14) 



We can rewrite much of this in terms of flexion. We define the 1- 
flexion as F = Fi + iF2, and the 3-flexion G = Gi -\- iG2, where 



F = -aaa>, 



G = -ddd^. 



(15) 
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F is manifestly spin 1 and G is spin 3. Comparing with equation 
JS]) and disentangling the individual components we find 

Fi = diti = di^i + 82^2, 

F2 — ^2/^ = ^172 — ^271, 

Gi = ^171 - ^272, 

G2 = ^172 + ^271- (16) 

Reorganising in terms of derivatives of shear, we can write 



^171 = 2^^' + ^') 
o>272 = \{Fi-Gi) 



^172 = -(^2 + ^2) 



^271 



-F2 + G2). 



(17) 



Hence we can write the D tensor in terms of the Pauli matrices as 
A,i = -Fil -\{Fi^ Gi)(j^ -\{F2 + G2)(Ji + \Cii(j2 

D^j2 = -F2I + ^(^2 - G2)CT3 - ^(^1 - G'Oai + ^C2iCT2 (18) 



where we have defined the turn, 

C^ = 2^^p 



(19) 



which in the complex notation we can write as C = 2dp, with C = 
Ci + iC2. This is a new distortion mode, which simply describes 
how the amount of image rotation in the Jacobian varies across the 
object. Like the rotation, it is not expected to be activated by gravity 
at our level of approximation. 

Separating into individual distortion components, we have 



-2D 



^jl5 



-2D 



.,2 J = Fi[2I + a3,ai]+F2[ai,2I -as] 
-\-Gi[a3, -cri] + G2[cri,a3] 
+CiHa2,0]+C2[0,-ia2]. (20) 

F, G and C provide six parameters for D. However, D has eight 
components, so it might initially be thought that there are eight de- 
grees of freedom in lensing distortions at this order. What do the 
remaining two parameters represent? We note that whereas C pre- 
multiplies cr2, F and G premultiply mixtures of /, ai and as. We 
can therefore seek a further mixture of these latter quantities. We 
can find this by writing the six known objects as 1-D lists of el- 
ements, treating these as vectors and seeking two further vectors 
which are orthogonal to these six and each other. Gaussian elimi- 
nation leads to the components 

[-2D^Jl, -2D^j2] = ... + Ti[-I + cTi + as, -I + ai - as] 

+T2[-I -ai+as,I + ai + as]. (21) 

where we have introduced the twist, T, which might appear to be 
another non-gravitational distortion mode. 

However, we note from equation iT3l that a lensed object has 
surface brightness at position Oi found using a second order term 
(A12 + Di2i)0i02/2; ie A12 and A21 do not occur indepen- 
dently of one another for any observational consequence. This sym- 
metrisation means that there are six, rather than eight, observational 
quantities at second order. We can cause a particular distorted im- 
age by either applying a twist or a turn. We will find the relation 
between the two below. 



We now have a complete list of distortions to second order. 
Explicitly, the full image distortion to this order is described by the 
A matrix. 



1-K 



+ 



1-K 

and the D tensor, 
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Gi G2 
G2 —Gi 

T1-T2 

Ti - T2 -2Ti - 2T2 

G2 —Gi 
-Gi —G2 

-2Ti + 2T2 Ti + T2 



T1+T2 







(23) 



From this list of the elements of D we can easily find how to con- 
vert between a twist mode and the turn which causes the same ob- 
servational consequences. If we have description a of a distortion 
with zero turn and non-zero twist {Ti,T2), this is equivalent to 
description b with zero twist, and turn given by 






1 -1 
-1 -1 



2 



(24) 



Equally if we start with description h, we can find description a 
using 

1/1 -I \ f Gl 






-1 



C? 



(25) 



The surface brightness mapping can be approximated to second or- 
der as 



fim = fs ( S^jOj + {A^j - d.j)Oj + -D^jkOjOk 



c^ fs{0^) + (A^j - 6^J)0J^^fs{0^) + -D,jkOjOkd^fs{0^). (26) 

Figure [1] gives an overview of the lensing quantities, their relation- 
ship to the gravitational potential and their transformation proper- 
ties. G is found by taking the derivative of p, but p and therefore 
G cannot be derived from the potential by taking derivatives. Thus 
p and C or T constitute additional degrees of freedom in the lens 
mapping beyond gravitational effects. 



4 BEHAVIOUR OF TWIST AND TURN 

4.1 Rotational Properties 

We can now examine how the new constituents of D transform 
under rotations R, by rotating a source-plane coordinate ft in the 
absence of A: 

Pi — RilDimnOmOn = RuDimnRmj ^jp^pRnk^kgOq 

= RilDimnRmjRnk^j^k (27) 

where primes denote rotated quantitites. But also if D' is the rotated 
tensor then ft = D'^jf^OjO'f^, so 



J-^ijk — ^ilJ-^lmni^mj^nk- 

Since we can write the rotation by angle (j) as 

_ / cos (j) sin q 
- sin 6 cos c 



(28) 



(29) 
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Q. 
C/3 

tensor 3 
tensor 2 

vector 1 
scalar 



potential 



deflection 



Jacobian 




12 3 

powers of differential operator 

Figure 1. Chart summarising the derivation of the lensing quantities by suc- 
cessive appHcation of the differential operators, and their respective trans- 
formation properties. 



we can write the transformation for D (and its constituent objects) 
as 



[D'iji,D'ij2] = [RilDlmlRmj COS + RilDlm2Rmj Sm (/ 
-RilDlmlRmj sin + RiiDlm2Rmj COS ( 



(30) 



If we set F = C = T = C2 = so that D only contains non-zero 
turn Ci, we find that this Ci object transforms into the equivalent 
C2 object after a 7r/2 rotation, and returns to its initial form after 
a 27r rotation; so Ci and C2 are the components of a vector. This 
is to be expected, as turn is the gradient of the scalar rotation field, 
and so is naturally a vector. 

Since a twist can be written in terms of a turn, we should ex- 
pect that twist will also be a vector. Indeed, if we set F = C = 
C = T2 = 0, the resulting Ti object transforms into the equivalent 
T2 object after a 7r/2 rotation, and returns to its initial form after a 
27r rotation; so Ti and T2 are also the components of a vector. 



4.2 Real Space Behaviour 

Now that we have established the rotational properties of the dis- 
tortions, we would like to visualise what effect they have on real 
images. We can use the first line of the mapping equation ([26]) to- 
gether with the D tensor of equation (l23l l to observe the effect of 
the second order image distortions on a Gaussian circular or ellip- 
tical image. 

Figure 12] shows the effect of 1 -flexion and 3-flexion. Here we 
have operated on objects with Gaussian surface brightness 



/(x) 



Aexp 



{x- 



(y-ycf 



2a2 



2a2 



(31) 



where (Tx = (Jy = 0.5" for the circular source, and ax = 
0.11", ay = 0.5" for the elliptical source. This gives an elliptic- 
ity e={al- al)/[al + al) = 0.9. 

Notice that 1 -flexion and 3-flexion affect the shape of both 
circular and elliptical objects in the figure. We can compare this to 
figure [3] which shows the effect of turn on circular and elliptical 
objects. Notice that turn has no discernible effect on the circularly 
symmetric source; we will show later that it indeed has zero effect 
on such a source. On horizontally or vertically elliptical objects, 
a pure Ci or C2 gives an arc which by eye appears similar to the 
impact of flexion, but it is truly a different mode of curvature with 
a distinct estimator which we will find below. 







Figure 2. Effect of 1 -flexion and 3-flexion on circular and elliptical (e 



0.9) Gaussian sources with a^ 



0.5". Top panel: unlensed objects; 



middle panel: Fi =0.2 arcsec ^; bottom panel: Gi =0.7 arcsec 



Figure |4] shows the equivalent effect of twist. It might appear 
from this figure that twist has a different effect to turn, but a com- 
bination of Ci and C2 can achieve the same effect as Ti or T2 ; it 
is just that a pure Ci distorts the object in a different direction to 
a pure Ti, for instance. Again, twist appears not to affect circular 
objects, and we will show this to be the case below. Its impact here 
is to turn horizontally or vertically elliptical objects into aerofoil 
shapes. Note the way in which T components engage with the el- 
lipticities to make aerofoils oriented in different directions. Positive 
Ti or negative T2 operate to twist horizontal objects into upward 
curving objects, with the front of the aerofoil pointing in opposite 
senses; negative Ti or positive T2 operate on horizontal objects to 
make downward curving objects. On the other hand, positive Ti 
and positive T2 bend vertical objects to the left, while negative Ti 
and negative T2 bend them to the right. 



4.3 Shapelet Space Behaviour 

We can gain insight into the behaviour of these distortion modes 
by examining their action in shapelet s pace. We use th e po- 
lar shapelets of Bernstein & JarvisI (l2002h . iRefregieJ ([20031), and 
iMassey & Refregier (2005). As is described in the latter paper, po- 
lar shapelets can be described by their number of radial nodes n and 
azimuthal nodes m, providing a basis set |n, m) for 2D localised 
objects. The shapelets require a length scale (3 to be set, which is 
the standard deviation of the zeroth shapelet, a 2-D circular Gaus- 
sian. Then an image |/) is the sum of the shapelets with appropriate 
coefficients: 



i/> = E/" 



(32) 
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Figure 3. Effect of turn on horizontally and vertically elliptical objects. Top 
panel: Ci = 1.3 arcsec"-*^ on a circular Gaussian with a = 0.5". Second 
panel: Ci = 1.3 arcsec"-*^ on elliptical Gaussians with (Tmajor = 0.5", 
ei = 0.9 and ei = —0.9. Third panel: Ci = —1.3 arcsec"^ on the 
same elliptical Gaussians. Fourth panel: C2 = 1.3 arcsec" ^ ; Bottom panel: 
C2 = — 1.3 arcsec"^. 



Figure 4. Effect of twist on horizontally and vertically elliptical objects. 
Top panel: Ti = 1.3 arcsec"^ on a circular Gaussian with a = 0.5". 



Second panel: Ti = 1.3 arcsec ^ on elliptical Gaussians an 



:0.5", 



ei = 0.9 and ei = —0.9. Third panel: Ti = —1.3 arcsec ^ on the 
same elliptical Gaussians. Fourth panel: T2 = 1.3 arcsec"^ ; Bottom panel: 
T2 = -1.3arcsec-^ 



and a lensed image is the result of applying various operators to the 
source: 

I/) = {l^Kk + pR + J^S^^F^F^ + G^G^+T^f^ + C^C^)\f){33) 

where the terms are for convergence, rotation, shear, 1 -flexion, 3- 
flexion, twist and turn respectively. We wish to discover what these 



operators are in terms of the ladder operators which act on the basis: 



.fi \ /n + m + 2, 
K\n,m) = y 1 



n + 1, ??i + 1) 



. I \ jn + m 

ar\n,m) = W — - — |n-l,m-l) 

-ti \ In — m-\-2 . 

a\\n,m) = W ^ |n + l,m- 

- I \ In — m . ^ ^. 

ai\n,m) = \ — - — |n-l,m + l). 



(34) 
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Figure 5. Effect of shear and flexion on a circular Gausian in shapelet space. 
Shapelet profiles are displayed for the real part of polar shapelets in the top 
part of the figure, and for the imaginary part in the bottom part of the figure. 
Shear (red) takes power from the Gaussian shapelet |0, 0} and places it in 
the spin-2 modes |2, ±2). 1-flexion (yellow) moves power to spin-1 modes, 
while 3-flexion (blue) moves power to the spin-3 modes |3, ±3). 



The ladder operators obey commutation relations 

ar,aj = 1 



\ai,ar 



\di,dl\ — \dl,dr\ — \dl,dl\ — 0. 



(35) 



We can write position and derivative operators in terms of these 
ladder operators: 



dx 

d_ 

dy 



l[dl 

IK 



+ dl -\- di -\- d^ 



■ d + dl 



ij -\- dl + dr 



a; + a + ai 



- ttr 



(36) 



Using these with equation (l26l ), we can find forms f or firs t and 
second order lensing operators. iMassev & Refregien (l2005h and 
iMassev et alJ (120070 have shown the forms for shear and flexion; 
these are summarised in figure \5\ which shows their effects on a 
circular Gaussian in shapelet space. Note that shear moves power 
from the |0, 0) mode to the spin-2 modes |2, ±2), while 1-flexion 
moves power to spin-1 modes, and 3-flexion moves power to the 
spin-3 modes |3, ±3). 

We can carry out similar calculations for twist and turn, using 
equation (l26l l together with equation (l36l ). For twist, after routine 



Figure 6. Effect of twist in shapelet space. Circular Gaussians are not af- 
fected, so here we show the power moved from the |2, 2} mode; note that 
twist pushes power into spin 1 and spin 3 modes. 



but extensive non-commutative algebra we find 



fi = -- M - djdi + dldrj dr{l - i) 
- M - dldr + djdi] di{l + i) 
—dl ( d\di — dldr ) (1 + 
—d\ { d]di — dldr ) (1 — i) 

T2 = -- il - d\di + dldrj dr{l +i) 

+ M - dldr + d\di\ di{l - i) 

+dl { d\di — dldr ) (1 — 

—d\ { d]di — dldr ) 0- + i, 



(37) 



where the factor of [3 takes into account the fact that the operators 
in equation (|36] ) work in units of (3. One can consider what happens 
to a circular (m = 0) source operated on by e.g. Ti; the third and 

fourth terms in the equation above for Ti vanish, as ( d\di —dldr] 
counts m; the first term initially acts with ar to move the state to a 
spin-1 state; then the term in brackets il — d\di + dldr ) operates 
to give zero. The second term similarly gives zero, resulting in twist 
having no effect on circularly symmetric objects. 

The effect of these operators is shown in figure [6] Since they 
have no impact for circular objects, we show the effect on the |2, 2) 
mode. Note that power is moved to neighbouring spin 1 and spin 
3 modes, with a rotation of (1 + i) or (1 — i) which gives the 
characteristic twisted form of the image. 

We carry out similar calculations for the equivalent turn de- 
scription, again using operators given by equation (|36] | together 
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Figure 7. Effect of turn in shapelet space. As for twist, circular Gaussians 
are not affected, so we again show the power moved from the |2, 2} mode. 
The Ci component (purple) pushes power into spin 1 and spin 3 modes on 
the opposite side of the diagram, while C2 (orange) does the same on the 
near side of the diagram. 



with the mapping equation (l26l ). We find 



Ci 



1(3 



f[-( 



1 + a\ai 



1 + alcir — ajcLl 1 ttr 



Idrj ai+ ( 
I ( didi — dldr ] -\- dl [ dldr — djdi J 



where f[i is the component after twist. Since the mean untwisted 
/ii is expected to be zero, we have the estimator 



T-est I -rriest 
1 -\-1-L2 



_2 i 77/4/11 



In a similar fashion, the turned /n coefficient is given by 

f[l=fll + ^{C\-iC2)f22 



(40) 



(41) 



so since the undistorted mean fu is expected to be zero, we obtain 
the estimator 



Cf * - iC| 



2V2i /i 
/3 /2: 



(42) 



Note the close relationship between the estimators for twist and 
turn, and the fact that they can indeed be written as superpositions 
of each other. 

However, it should be noted that these will only be pure es- 
timators for twist/turn if 1 -flexion is absent (or negligible). This 
is because 1 - flexion also moves power into the fu mode (see 
iMassev et al.! (l2007h ) in such a fashion that our estimator (e.g. for 
twist) is truly 



i7r/4fn_ 



(i^ 



/22 



" 2/22 



(3/00—3/40) 



F*e^ 



rr/4 



4/22 



-(y2/22-3y6/42)(43) 



The estimator is still of value to us despite this complication, as we 
wish to use it to see if there is a twist-like systematic in our survey; 
we now see that, at this shapelet order, systematic flexion can gen- 
erate a twist-like effect. This simple estimator is therefore a test of 
combined second-order systematics, or of (real or systematic) twist 
on scales where flexion is negligible. 

Nevertheless, it should be kept in mind that a pure estimation 
of twist will require a more extensive joint chi- squared fit of twist 
and flexion to several further orders of shapelets, in order to fully 
remove the degeneracy. 



C2 — —— \—[l + d\di— dldr \ dl — il + dldr — d\di j dr 



+ 



d\ I djdi — dldr ] + dl I dldr — djdi J (38) 



Here again we find that the impact on m = states is zero, using an 
identical argument to above. The effect of these operators is shown 
in figure |7] Again we show the effect on the |2, 2) mode; as for 
twist, power is moved to neighbouring spin 1 and spin 3 modes, 
but there are a different range of activated modes for Ci and C2, 
due to the different factors of i. 

Now that we can describe twist/turn in shapelet space, we are 
in a position to construct practical estimators for measuring these 
quantities. 



5 CONSTRAINING TWIST AND TURN 
5.1 Simple Estimators 

We can use the ladder operator form for twist to find a simple esti- 
mator for the new distortions. We consider the power that finishes 
in the fu component. 



/n=/ii-fe-^^/'(Ti+iT2)/22 



(39) 



5.2 Correction for Centroid Shift 

In addition, these simple estimators need cor recting for the fact 
that tw ist/turn moves the centroid of the object. iGoldberg & BaconI 
(l2005h show that the centroid is moved by the D tensor according 
to 



AOi 



^l) ( 1^111 + ^^212 + Id221 



9162) {D112 + D12I + i^222) - {02)-Di22 



AO2 = -{Ot)-D2ii - {O1O2) {D221 + D212 + Dm) 

-{Ol) (\d222 + \Di21 + \d^12 



(44) 



where we have written a form which assumes less symmetry than 
Goldberg & Bacon; this is necessary for our generalised D tensor. 
Putting the values of the D tensor, equation (l23l l, into this equation 
we find 

o2 

A6>i+iA6>2 ^ Jn[^F + 5F*e + Ge* + iC*e + (i - l)Te](45) 

where R^ i s the size quad rupole and e is the unweighted ellipticity 
as given in lMassev et al.l ([2007). Note again that for circular (e = 
0) objects, twist and turn have no effect. 
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iMassev et alJ ("2007) showed that the effect of the centroid 
shift is to alter observable flexion estimators Fi + iF2 by substract- 
ing a term {AODr + AO* Di), and the same applies here. The shift 
operators D are given by 



Dr = - (al - an ^^ ^ 9 (^^^ ~ ^"^ ) 

Hence we find that equation ( [39l becomes 

111 = 111 - 2^ ^/22 

- -^T [{i - l)(ei + ie2){foo - /20) 
+ V2{i + l){ei-ie2)f22\ 



(46) 



(47) 



By dividing both sides by /22, we can therefore propose the cor- 
rected estimator 

^est ^ _lg^^/4/il A _ / R^jei +ie2)(/oo -/2o) \\ .^g. 



If we label the term in brackets as B, we similarly find for the turns 



C 



2V2i /i 



(49) 



B(3 /22 

So our corrected estimators differ from our naive estimators only 
by a factor of B. Again, twist can be written as a superposition 
of turns, and vice versa. As in the previous section, the presence 
of flexion would mean that these are not pure estimators of twist; 
the estimators should again be interpreted as showing combined 
second-order systematics, or (real or systematic) twist on scales 
where flexion is negligible. 



6 FIRST MEASUREMENTS WITH TWIST/TURN 
ESTIMATORS 

We are now in a position to measure these twist/turn estimators on 
real data. We use the STAGES niosaic observed with th e Hubble 
Space Telescope JGrav et al.ll2007l : iHevmans et al.ll2008h . This is a 
0.25 square degree field observed with the Advanced Camera for 
Surveys in the F606W band, covering 80 ACS tiles in 80 orbits. 
Drizzling is used to obtain an effective pix el size of 0.03" . 

We use the same galaxy catalogue as iHevmans et al.l (l2008h . 
deconvolving and de composing all objects into shapelets using th e 
methods developed in'RefregieiJ (l2003h : lRefregier & Baconl(l2003h : 
iMas sev & Refregier (2005). The analysis will be described in full 
in Bacon et al (2009); we obtain a shapelet catalogue for 56,000 
galaxies, together with measures of f3 and R^ for all objects. The 
shapelets are normalised so that /oo = 1- 

We measure the twist/turn estimators for objects with F606 
magnitude < 23.5 using equations (|48] | to (|49] |; we find that for 
this sample, B — — 1.59 =b 0.01. Since the twist and turn measure- 
ments are equivalent, here we choose to present results in terms of 
twist. The histogram of twist estimators is shown in figure [H this 
includes 3cr cuts for outliers with \T\ > 6.0 arcsec"^, and we only 
consider objects with f3 >\ pixel to avoid oversampling. The first 
thing to note is that our twist estimator is more noisy than shear 
and flexion estimators, having a standard deviation in one compo- 
nent of 2.0 arcsec"^; much of this noise is due to intrinsic shape 
variance of the objects. The turn estimator as defined has a larger 
standard deviation of 2.9 arcsec"^. 

We find mean values over the STAGES survey of Ti = 




-6-4-2 2 4 6 

T (arcsec"^) 

Figure 8. Histogram of twist estimator values in one component in the 
STAGES survey, for objects with F606 magnitude< 23.5, \T\ < 6.0 
arcsec"^ and f3 > 1 pixel. 



-0.016±0.036arcsec"\T2 = -0.009 ± 0.037 arcsec"\ These 
are consistent with zero, as we might hope for a systematic mode. 
At present the constraint is fairly weak, as gravitational flexion sig- 
nals are at the level of 0.001 to 0.01 arcsec"^; however, in up- 
coming lensing surveys the much larger area will lead to twist/turn 
constraints at the 10~^ level, which will provide important checks 
on systematics. 

We can further explore whether the twist/turn estimator is ac- 
tivated as a systematic in the STAGES survey by measuring its cor- 
relation functions. As with shear correlation functions, the twists 
should be rotated before they are correlated; however, while shear 
has to be rotated by factors of e*^'^ where is the position angle of 
the line joining a pair of objects, twist has to be rotated by factors 
of e*"^ on account of its vector nature: 



T[°* = Ti cos + T2 sin 
T2°* = -Ti sin + T2 cos 

We can then construct correlation functions 

cliio) = {Tr\e,)Tr\e, + 0)) 



(50) 



cUo) 
cUo) 



{n''\e^)n''\e. + e)) 



in 



)n''\o, + o)) 



(51) 



We have measured these correlation functions for twist estimators 
in STAGES, and display the results in figure |9] Here error bars 
are estimated by cr^/y^Np^ where a is the standard deviation 
of twist and A^pairs is the number of galaxy pairs in a bin. We find 
that the correlation functions are almost all consistent with zero 
signal, with reduced xl = 0.87, 0.53, 0.39 for Cfi, CJ2, ^^2 re- 
spectively. 



7 CONCLUSION 

In this paper, we have written down for the first time the full weak 
image distortion relevant to weak lensing, to second order. This 
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Figure 9. STAGES twist estimator correlation function. Solid line: (Ti Ti } 
correlation function. Dashed line: {T2T2) correlation function. Dotted line: 
{T1T2} correlation function. 



involved the discovery of a new image distortion, the twist, which 
can be written in an alternative form as the turn. 

We reviewed weak lensing distortions to first order, recouch- 
ing the lens mapping in terms of Pauli matrices. We noted the ex- 
istence of a non-gravitational mode, the rotation. This sets a prece- 
dent which we also see at second order. 

We then extended the formalism to second order; at this point 
it became clear that the gradient of the rotation gives a new mode, 
which we call the turn. A further orthogonal mode in D was found 
by seeking a combination of Pauli matrices orthogonal to all known 
modes; this new mode was called the twist. We showed that twist 
and turn can be written in terms of each other regarding their ob- 
servational consequences, and were then able to write down the full 
image distortion mapping to second order. 

We explored the properties of twist/turn, finding that it is a 
vector quantity. Its visual effect was shown, as was its impact in 
shapelet space. With the ladder operator formalism we showed that 
twist/turn has no effect on circularly symmetric objects, but only 
objects with non-zero ellipticity. We saw how twist/turn moves 
power from spin-2 modes to modes with spin-1 and spin-3. 

Using our ladder operator forms for twist and turn, we found 
simple estimators for the distortions; however, twist/turn causes a 
centroid shift which needs to be taken into account. This leads to 
the inclusion of a common factor B in the estimators for twist and 
turn. In addition, any flexion present would contribute to our simple 
estimator, which should therefore be used either as a means to de- 
tect any second-order systematic, or to measure real or systematic 
twist on scales where flexion is negligible. 

We used these estimators to constrain twist for the first time 
in the HST STAGES survey. We noted that our estimator has a 
larger intrinsic noise scatter compared to 1- or 3-flexion, but its 
mean value across the survey is already at an interesting level for 
checking large flexion systematics. We found that mean twist is 
consistent with zero in this survey. We also measured twist correla- 
tion functions, and found that they too were consistent with zero. 

The two quantities introduced in this study complete the set 



of distortions to second order. They will be of use in testing for 
systematic effects, and have a certain elegance of their own. 
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